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Abstract 



In this paper we study different algorithms for reflected backward stochastic differential equa- 
tions (BSDE in short) with two continuous barriers based on binomial tree framework. We intro- 
duce numerical algorithms by penalization method and reflected method respectively. In the end 
simulation results are also presented. 

Keywords: Backward Stochastic Differential Equations with two continuous barriers, Penal- 
ization method, Discrete Brownian motion, Numerical simulation 

AMS: 60H10, 34K28 

1 Introduction 

Non-linear backward stochastic differential equations (BSDEs in short) were firstly intro- 
duced by Pardoux and Peng ([21J, 1990), who proved the existence and uniqueness of the 
adapted solution, under smooth square integrability assumptions on the coefficient and the 
terminal condition, plus that the coefficient g(t, u, y, z) is (t, a;)-uniformly Lipschitz in (y, z). 
Then El Karoui, Kapoudjian, Pardoux, Peng and Quenez introduced the notion of reflected 
BSDE (RBSDE in short) QllJ, 1997) with one continuous lower barrier. More precisely, a 
solution for such an equation associated to a coefficient g, a terminal value £, a continuous 
barrier L t , is a triplet (Y t , Z t , A t ) < t <T of adapted processes valued in R 1+d+1 , which satisfies 
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and Y t > L t a.s. for any < t < T. A t is non-decreasing continuous, and B t is a d- 
dimensional Brownian motion. The role of A t is to push upward the process Y in a minimal 
way, in order to keep it above L. In this sense it satisfies f Q (Y s — L s )dA s = 0. 

Following this paper, Cvitanic and Karatzas ([9], 1996) introduced the notion of reflected 
BSDE with two continuous barriers. In this case a solution of such an equation associated 
to a coefficient g, a terminal value £, a continuous lower barrier L t and a continuous upper 
barrier U t , with L t < U t and L T < £ < Ut a.s., is a quadruple (Y t , Z t , A t , K t ) < t < T of adapted 
processes, valued in R 1+d+1 , which satisfies 



and L t < Y t < U t , a.s. for any < t < T. Here A t and K t are increasing continuous process, 
whose roles are to keep the process Y between L and U in such a way that 



In view to prove the existence and uniqueness of a solution, the method is based on a 
Picard-type iteration procedure, which requires at each step the solution of a Dynkin game 
problem. Furthermore, the authors proved the existence result by penalization method when 
the coefficient g does not depend on z. In 2004 Q16J), Lepeltier and San Martin relaxed in 
some sense the condition on the barriers, proved by a penalization method an existence result, 
without any assumption other than the square integrability one on L and U, but only when 
there exists a continuous semi-martingale with terminal value £, between L and U. More 
recently, Lepeltier and Xu ([IB]) studied the case when the barriers are right continuous 
and left limit (RCLL in short), and proved the existence and uniqueness of solution in 
both Picard iteration and penalization method. In 2005, Peng and Xu considered the most 
general case when barriers are just L 2 -processes by penalization method, and studied a 
special penalization BSDE, which penalized with two barriers at the same time, and proved 
that the solutions of these equations converge to the solution of reflected BSDE. 

The calculation and simulation of BSDEs is essentially different from those of SDEs (see 
[T4]). When g is linear in y and z, we may solve the solution of BSDE by considering its 
dual equation, which is a forward SDE. However for nonlinear case of g, we can not find the 
solution explicitly. Here our numerical algorithms is based on approximate Brownian motion 
by random walk. This method is first considered by Peng and Xu [25]. The convergence 
of this type of numerical algorithms is proved by Briand, Delyon and Memin in 2000 ([!]) 
and 2002 [5] . In 2002, Memin, Peng and Xu studied the algorithms for reflected BSDE 
with one barrier and proved its convergence (cf. [20J ) . Recently Chassagneux also studied 
discrete-time approximation of doubly reflected BSDE in [6]. 

In this paper, we consider different numerical algorithms for reflected BSDE with two 
continuous barriers. The basic idea is to approximate a Brownian motion by random walks 
based on binary tree model. Compare with the one barrier case (cf. [2U]), the additive barrier 
brings more difficulties in proving the convergence of algorithm, which requires us to get finer 
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estimation. When the Brownian motion is 1- dimensional, our algorithms have advantages in 
computer programming. In fact we developed a software package based on these algorithms 
for BSDE with two barriers. Furthermore it also contains programs for classical BSDEs and 
reflected BSDEs with one barrier. One significant advantage of this package is that the users 
have a very convenient user-machine interface. Any user who knows the basics of BSDE can 
run this package without difficulty. 

This paper is organized as follows. In Section 2, we recall some classical results of reflected 
BSDE with two continuous barriers, and discretization for reflected BSDE. In Section 3, we 
introduce implicit and implicit-explicit penalization schemes and prove their convergence. 
In Section 4, we study implicit and explicit reflected schemes, and get their convergence. In 
Section 5, we present some simulations for reflected BSDE with two barriers. The proof of 
convergence of penalization solution is in Appendix. 

We should point out that recently there have been many different algorithms for com- 
puting solutions of BSDEs and the related results in numerical analysis, for example [3], [1], 
[Z] 5 [E], [I0], [13], [19], [26]. In contrast to these results, our methods can easily be realized 
by computer in 1-dimensional case. In the mult i- dimensional case, the algorithms are still 
suitable, however to realize them by computer is difficult, since it will require larger amount 
of calculation than 1-dimensional case. 

2 Preliminaries: Reflected BSDEs with two barriers 
and Basic discretization 

Let (fljJ 7 , P) be a complete probability space, (B t )t>o a 1-dimensional Brownian motion 
defined on a fixed interval [0, T], with a fixed T > 0. We denote by {^t}o<t<T the natural 
filtration generated by the Brownian motion B, i.e., Tt = cr{B s ;0 < s < t} augmented 
with all P-null sets of T. Here we mainly consider 1-dimensional case, since the solution of 
reflected BSDE is 1-dimensional. In fact, we can also generalize algorithms in this paper to 
multi-dimensional Brownian motion, which will require a huge amount of calculation. We 
introduce the following spaces for p e [1, oo): 

• L p (jF t ) :={R-valued jF t -measurable random variables X s. t. < oo}; 

• L^-(0, t) :={M- valued and jF t -adapted processes ip defined on [0, t], s. t. E J Q * \Lp s \ p ds < 
oo}; 

• S p (0,£) :={IR-valued and jF t -adapted continuous processes tp defined on [0,£], s. t. 
£[sup <t< T \<ft\ 2 ] < oo}; 

• A p (0,t) :={increasing processes in S p (0,t) with A(0) = 0}. 
We are especially interested in the case p = 2. 
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2.1 Reflected BSDE: Definition and convergence results 

The random variable £ is considered as terminal value, satisfying £ G L 2 (jF r ). Let g : 
[0, T] x R x R — ► R be a ^-uniformly Lipschitz function in (y,z), i.e., there exists a fixed 
/i > such that 

\g(t,y 1 ,zi)-g(t,y 2 ,z 2 )\ < fi(\yi - y 2 \ + \zi - z 2 \) (1) 

Vt G [0,T],V(yi,2:i), (7/2,^2) G R x R. 

And gr(-, 0, 0) is square integrable. 

The solution of our BSDE with two barriers is reflected between a lower barrier L and 
an upper barrier U, which are supposed to satisfy 

Assumption 2.1 L and U are Tt-progressively measurable continuous processes valued in 
R, such that 

E{ sup ((L t ) + ) 2 + sup ((t^n 2 ] < 00. (2) 

0<t<T 0<t<T 

and there exists a continuous process X t = X — f Q a s dB s + V t + — Vf where a G L^-(0, T), 
V + and V~ are (Ft) -adapted continuous increasing processes with E\\V^\ 2 ] + E\\V^\ 2 ] < 00 
such that 

L t <X t < U t , P-a.s. forO<t<T. 

Remark 2.1 Condition (0|) permits us to treat situations when Ut = +00 or L t = — oo ; 
t G [0,T], in such cases the corresponding reflected BSDE with two barriers becomes a 
reflected BSDE with a single lower barrier L or a single upper barrier U , respectively. 

Definition 2.1 The solution of a reflected BSDE with two continuous barriers is a quadruple 
(Y,Z,A,K) G S 2 (0,T) x 14.(0, T) xA 2 (0,T) x A 2 (0,T) defined on [0,T] satisfying the 
following equations 

-dY t = g(t,Y t , Z t )dt + dA t - dK t - Z t dB t , Y T = £ (3) 
L t < Y t < U t , dA t > 0, dK t > 0, dA t ■ dK t = 0. 

and the reflecting conditions 

[ {Yt-L t )dA t = [ (Y t -Ut)dK t = 0. (4) 
Jo Jo 

To prove the existence of the solution, penalization method is important. Thanks to the 
convergence results of penalization solution in [16], [15] for continuous barriers' case and 
methods in [23j, we have the following results, especially it gives the convergence speed of 
penalization solutions. 
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Theorem 2.1 (a) There exists a unique solution (Y, Z, A, K) of reflected BSDE, i.e. it 
satisfies (dj), Moreover it is the limit of penalization solutions (Y™' 10 , Z™ ,p , A™' p , K™' p ) 
as m — > oo then p — > oo, or equivalent as q — > oo then m — > oo. Here the penalization 
solution (l^ m ' p , Z™ ,p , A™' p , K™' p ) with respect to two barriers L and U is defined, for m G N, 
pGN, as the solution of a classical BSDE 

-d? t m ' p = g(t,Y t m ' p ,Zr ,p )dt + m(Y t m ' p - L t y dt - q(Y t m > p - U t ) + dt - Z™' p dB tl (5) 
Y^ p = £. 

And we set A™' p = m j*{Y s m ' p - L s )'ds, K™' p = p J*(Y s m 'P - U s ) + ds. 

(b) Consider a special penalized BSDE for the reflected BSDE with two barriers: for any 
p E N, 

-dY p = g{t,Y p ,Z p )dt + p{Y p - L t ydt-p(Y p -U t ) + dt- Z p dB t , (6) 

with A P = J ' p(Y p — L s )~ds and K p = J Q * p(Y p — U s ) + ds. Then we have, as p — > oo, Y t p — > Y t 
in S 2 (0,T) ; Z p -> Z t in l£(0,T) and A p -> A t weakly in S 2 (0,T) as well as K p -> 
Moreover there exists a constant C depending on £, ^(i, 0,0), /i, L and J7, sttc/i that 



E[ sup |y/ - Y t \ 2 + / T (|Zf - Z t | 2 )dt + sup [(A t - K t ) - (A p - K p )f] < SL. 

0<t<T Jo 0<t<T y/P 

The proof is based on the results in [16] and [23] , we put it in Appendix. 



(7) 



Remark 2.2 In the following, we focus on the penalized BSDE as (0), which consider the 
penalization with respect to the two barriers at the same time. And p in superscribe always 
stands for the penalization parameter. 

Now we consider a special case: Assume that 
Assumption 2.2 L and U are ltd processes of the following form 

L t = L + [ l s ds+ [ a l s dB s , (8) 
Jo Jo 

U t = U + u s ds+ / a u s dB s . 
Jo Jo 

Suppose that l s and u s are right continuous with left limits (RCLL in short) processes, a[ 

p T 2 / 1 2 2 2 

and a" are predictable with E[J Q [\l s \ + \a s \ + \u s \ + |cr"| )ds < oo. 

It is easy to check that if L t < U t , then Assumption 12.11 is satisfied. We may just set 
X = L or U, with a s = o~\ or cr" and V ± = L ifds or J Q ufds. Here If (resp. uf) is the 
positive or the negative part of I (resp. u). As Proposition 4.2 in [llj, we have following 
proposition for two increasing processes, which can give us the integrability of the increasing 
processes by barriers. 
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Proposition 2.1 Let (Y, Z, A, K) be a solution of reflected BSDE Then Z t = a\, a.s.- 

dP x dt on the set {Y t = L t }, Z t = o~f, a.s.-dP x dt on the set {Y t = U t }. And 

< dA t < l {Yt =L t }[g(t,Lt,a l t ) + l t ]~dt, 
< dK t < l {Yt=Ut} [g(t,U t ,o-?)+u t ] + dt. 

So there exist positive processes a and [3, with < a t , f3 t < 1, such that dA t = atl{Y t =L t }[g(t, L t , o~ l t ) + 
l t ]-dt, dK t = p t l {Yt=Ut} [g(t, U t , <r«) + u t }+dt. 

Proof. We can prove these results easily by using similar techniques as in Proposition 
4.2 in [11], in view that on the set {L t = Ut}, we have a\ = cr" and l t = Ut- So we omit the 
details of the proof here. □ 

In the following, we will work under Assumption 12.21 

2.2 Approximation of Brownian motion and barriers 

We use random walk to approximate the Brownian motion. Consider for each j = 1, 2, • • • , 

WS] 

B n t := 41 V e?, for all < t < T, 5 = -, 

J n 

where {e"}™ =1 is a {1, — l}-valued i.i.d. sequence with P(Bj = 1) = P (e" = —1) = 0.5, i.e., 
it is a Bernoulli sequence. We set the discrete filtration := a{e 1 , ■ ■ ■ , e"} and tj = jS, for 
< j < n. We denote by T) t the space of RCLL functions from [0,t] in IR, endowed with 
the topology of uniform convergence, and we assume that: 

Assumption 2.3 T : D^-> R is K-Lipschitz. We consider £ = T(B), which is Tt- 
measurable and £ n = T(B n ), which is Q™-measurable, such that 

£[|£| 2 ]+sup£[|ff] <oo 

n 

Now we consider the approximation of the barriers L and U. Notice that L and U are 
progressively measurable with respect to the filtration (J 7 *), which is generated by Brownian 
motion. So they can be presented as a functional of Brownian motion, i.e. for each t G [0, T], 
L t = {B s )o< s <t) and U t = $ 2 {t, (B s ) < s < t ), where •) and * 2 (*, •) : D*^ R. And 
we assume that ^i(t, •) and ^2(^,-) are Lipschitz. Then we get the discretizaton of the 
barriers L] = * 1 {t j , {B^) < s < t ) and U? = * 2 (^, (5 s w ) < s < t ). If L t < U t , then L j < U r On 
the other hand, we mainly consider barriers which are Ito processes and satisfy Assumption 
12.21 So we have a natural approximation: for 1 < j < n, 

j-i 3-1 
L] = L + Sj2h + J2 al i £ i+i^^ 

i-i i-i 
U» = Uo + S^ + zZ^ 6 ^ 6 

i=0 i=0 
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where k = l ti , o\ = a l t ., Ui = u ti , erf = af.. Then L™ and U? are discrete versions of L and 
U, with supJ J E[sup j ((L™)+) 2 + sup 3 -((£/™)-) 2 ] < oo and L" < L/t* still hold. In the following, 
we may use both approximations. 

In this paper, we study two different types of numerical schemes. The first one is based 
on the penalization approach, whereas the second is to obtain the solution Y by reflecting 
it between L and U and get two reflecting processes A and K directly. Throughout this 
paper, n always stands for the discretization of the time interval. And process (0™ )o<j< n is 
a discrete process with n + 1 values, for <fi = L, U, y p , y, etc. 



3 Algorithms based on Penalization BSDE and their 
Convergence 

3.1 Discretization of Penalization BSDE and Penalization schemes 

First we consider the discretization of penalized BSDE with respect to two discrete barriers 
L n and U n . After the discretization of time interval, we get the following discrete backward 
equation on the same interval [tj, tj+i], for < j < n — 1, 

VT = y^x + gi^yf^z^S + a^-k^-z^VSe]^, (9) 
af n = P 5(y?> n -L])-, ^ n = P 6{yf n -U-y. 

The terminal condition is y p ' n = £ n . Since for a large fixed p > 0, (EJ) is in fact a clas- 
sical BSDE. By numerical algorithms for BSDEs (cf. [21]), explicit scheme gives Zj' n = 
27s^y?+'^ £ i= 1 ~ £/j+iUj=— 1)3 an d y P j ,n is solved from the inversion of the following mapping 

where Q p (y) = y - g(t v y, zf n )6 - p5(y - L])~ + P 5(y - U?) + , 

by substituting E{y v £ x \Q™\ = §(yJ£| E » +1 =i + yf^\ e" +1 =-i) into h. And increasing processes 
a P ' n and H' n will be obtained from ([9j). 

In many cases, the inversion of the operator O p is not easy to solve. So we apply the 
implicit -explicit penalization scheme to ([HD, replacing y p ' n in g by and get 



a' 



p,n r/— p,n r"\ — T"Pi n c/—p,n 



J 



In the same way, we get zf n = ^(Vj+i\^=i ~ Vj+iU^-i) an d 

vT = m^m + 9(t jt E[f£mur) 5 + " ! r - v 'i"- no.) 

Solving this equation, we obtain 

vT = m^m+g{^E[^ x \g%z^)5 



with E]^'^\Qf\ = |(^+i| e y +1 =i + Vj+i\e^ +1 =-i)- For increasing processes, we can get them 
from 

3.2 Convergence of penalization schemes and estimations 

First we give the following lemma, which is proved in [20]. This Gronwall type lemma is 
classical but here it is given with more detailed formulation. 

Lemma 3.1 Let a, b and a be positive constants, 5b < 1 and a sequence (i>j)j=i ) ...n of positive 
numbers such that, for every j 



v j 



+ a < a + bS^Vi. (11) 



Then 

sup Vj + a < a£s(b), 

where £&{b) = 1 + Y^pLx (1 + ■ • ■ (1 + (p — l)S), which is a convergent series. 

Notice the £s{b) is increasing in 5 and 5 < |, so we can replace the right hand side of 
( TTTj) by a constant depending on b. 

We define the discrete solutions, (Y p ' n , Z p,n , A P ' n , Kf' n ) by the implicit penalization 
scheme 

[t/S] [t/S] 
VP, n _ „.P,n *7P,n _ p,n A p,n _ p,n jfP,n _ ip,n 

X t — y\t/S\i — z [t/6}i *H — ^% ) — rn ) 

m=0 m=0 

or (Yf' n , Zf' n , ~A*t' n , K p ' n ) by the implicit-explicit penalization scheme, 

_ WS] [t/s] 

yP>» _ tf>,n yP,n _ -p,n ~TP, n _ — p, n zfP,n _ TP< n 

I t — y[t/s\> *t — z [t/sY * ~ 2.^1 m ' * ~~ m ' 

m=0 m=0 

Let us notice that the laws of the solutions (Y p , Z p , A p , K p ) and (Y p > n , Z p > n , A p ' n , K p > n ) or 
(Y p ' n , Z p ' n ,T' n , K p ' n ) to penalized BSDE depend only on (P Bj r-\P B ),9, ^i 1 (Pb)^ 2 1 ( p b)) 
and (Pen, r _1 (P B «), g, (Pfin), ^2 1 (Pb™)) where P B (resp. Pb™) is the probability intro- 
duced by B(resp. and / _1 (P.b) (resp. /"'(Pb^)) is the law of f(B) (resp. f{B n )) for 



S 



/ = r, ^i, So if we concern the convergence in law, we can consider these equations on 
any probability space. 

By Donsker's theorem and Skorokhod representation theorem, there exists a probability 
space, such that sup 0<t<T \B™ — B t \ — ► 0, as n — ► oo, in L 2 (jF T ), since Sk is in L 2+<5 . So we 
will work on this space with respect to the filtration generated by B n and B, trying to prove 
the convergence of solutions. Thanks to the convergence of B n , (L n ,U n ) also converges 
to (L,U). Then we have the following result, which is based on the convergence results 
of numerical solutions for BSDE (cf. [I], [5]) and penalization method for reflected BSDE 
(Theorem 12.11) . 

Proposition 3.1 Assume \2.3\ holds, the sequence (Y p ' n ,Zf' n ) converges to (Y t , Z t ) in fol- 
lowing sense 

r T 

lim lim E[ sup \Y p ' n - Y t \ 2 + / \Z p ' n - Z s \ 2 ds] -> 0, (12) 

p-»oon-*oo <t<T Jo 

and for < t < T, A P ' n — K p ' n — > A t — K t in L 2 (jF t ) ; as n — > oo, p — > oo. 
Proof. Notice 

E[sup \Y t p > n -Y t \ 2 + / \Z P / 1 - Z s \ 2 ds] < 2£[sup \Y t p ' n - Y t p \ 2 + / \Z p ' n - Z p \ 2 ds] 

0<t<T Jo 0<t<T Jo 

+2£[sup \Y P -Y t \ 2 + / \Z P -Z s \ 2 ds}. 

0<t<T Jo 

By the convergence results of numerical solutions for BSDE (cf. [I], [5]), the first part tends 
to 0. For the second part, it is a direct application of Theorem 12.11 of the penalization 
method. So we get (1121) . For the increasing processes, we have 

E[((A p t ' n - K p ' n ) - (A t - K t )) 2 } < 2E[((Ar - KD - (A p - K p )) 2 ] 

+2E[((A p - K p ) - (A t - K t )) 2 ] 

< 2E[((Ar - KD - (A p t - K p )) 2 ] + 

in view of (JTj) . While for fixed p, 

A p ' n -Kr = Y£' n -Yf' n - ^ g(s,Y p>n ,Z p ' n )ds+ f Z p s ' n dB™, 

Jo Jo 

A p -K p = Y p -Y p - f g(s,Y p ,Z p )ds+ f Z p dB s , 

Jo Jo 

from Corollary 14 in [5], we know that J ZP' n dB™ converges to j Q Z^dB s in S 2 (0,T), as 
n — > oo, then with the Lipschitz condition of g and the convergence of Y p,n , we get (A P,n — 
K p ' n ) -> (A t - K t ) in L 2 (JS), as n -> oo, p -> oo. □ 

Now we consider the implicit-explicit penalization scheme. From Proposition 5 in [25] , we 
know that for implicit-explicit scheme, the difference between this solution and the totally 
implicit one depends on \i + p for fixed p e N. So we have 
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Proposition 3.2 For any peN, when n — > oo, 

E[sup |T^ ,n -y t p '"| 2 + f \z^ n - Zf"| 2 ds^O, 

0<t<T Jo 

with (A 9 ^ - 7Yf n ) - (A p t ' n - K p ' n ) ^Oin L 2 (F t ), forO<t<T. 

Proof. The convergence of (Y t ' ,~Z%' n ) to (Y p,n ,Z p,n ) is a direct consequence of Propo- 
sition 5 in [23]. More precisely, there exists a constant C which depends only on /x and T, 
such that 

£[ sup IFf " - Yf n | 2 ] + £ / |^' n - Z p ' n \ 2 < C<5 2 . 

0<t<T JO 

Then we consider the convergence of the increasing processes, notice that for < t < T, 

9 (s,Y P ; n X' n )ds+ fz^oIB:, 
Jo 

compare with Af n - Kf n = Y p ' n - Y t p ' n - / Q * g(s, Yf> n , Z^ n )ds + f Q Z^ n dB n s , thanks to the 
Lipschitz condition of g and the convergence of (Y ,Z ), we get A? -iff - ,4? n -lf? n , 
in L 2 (jF t ), as n — > oo, for fixed p. So the result follows. □ 

Remark 3.1 From i/ws proposition and Proposition \3.1l we get the convergence of the 
implicit-explicit penalization scheme. 

Before going further, we prove an a-priori estimation of (y p,n , z p ' n , a p,n , k p,n ). This result 
will help us to get the convergence of reflected scheme, which will be discussed in the next 
section. Throughout this paper, we use C^... to denote a constant which depends on (p, if), 
■ ■ ■. Here 0, ip, ' ' ' can be random variables or stochastic processes. 

Lemma 3.2 For each p e N and 5 such that 6(1 + 2/j, + 2/i 2 ) < 1, there exists a constant c 
such that 

n 1 n 1 n 

£[-p m 2 + E l*rf * + ^E K'f + - 6 E l*ff ] < c c^^- 

J 3=0 ^ 3=0 P ° 3=0 

idere C^n^^n^n depends on £ n ; g(t, 0,0), (-^ n ) + and (U n )~ , while c depends only on \x and 
T. 

Proof. Recall (jUJ), we apply 'discrete Ito formula' (cf. [20]) for (y p ' n ) 2 , we get 

71— 1 to—I 

£[|yTf+Ei*n 2 *] < ^[iei 2 +2Eyrb(*^r^ni^] 

i=3 i=j 

n-1 

i=j 



n - n = Yl' 11 -T t ' n - ( 

Jo 
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Since yf' n • af n = -p5{{yf n - £™)~) 2 + pttW" - L™)~ = + L>f » and »f B ■ Af » = 

PW - KD + ) a + U?p5( V r - Un + = js k ? n + W". we have 

1 n—1 1 n—1 1 n—1 

E U n \ 2 + 2 E i^i 2 5] + 2^[- $>n 2 + E(*n 2 ] 

n—1 n—1 n—1 n—1 

< E[\C\ 2 + E °> 0)| 2 5 + (1 + 2^ + 2^ 2 ) £ |yn 2 5 + 2 »;V»r + 2 53(&T)-Af "] 

i=j i=j i=j i=j 

n—1 n—1 n—1 

< e[\c\ 2 + E fe> °> °)i 2 ^] + (! + + 2/i 2 )^ E itf'T + -^(E a H 2 

i=j i=j i=j 

n—1 

sup m) + ) 2 + -E(J2kD 2 + snp ((^) + ) 2 - 

j<i<n— 1 P • • j<i<n— 1 
«=j 

Since L n and t/ n are approximations of Ito processes, we can find a process X™ of the form 
X™ = Xq — J2i=o a i £ ?+i VS+V J ^~ n — Vf n , where V^ n are (/"-adapted increasing processes with 
-^[IK + ™| 2 + lKT n | 2 ] < +°°) an< i Ej ^ X™ < U" holds. Then we apply similar techniques of 
stopping times as in the proof of Lemma 2 in [16J for the discrete case with L™ < X™ < 
we can prove 

n—1 n—1 n—1 

£(X>n 2 + ^(E^ ,n ) 2 < Mc^, 9 , X n +Ej2[\ y r\ 2 + \# n \*\s). 

i=j i=j i=j 

While X n can be controlled by L n and U n , we can replace it by L n and U n . Set a = (3 = 12/j, 
in the previous inequality, with Lemma 13.11 we get 

n—1 1 n—1 1 n—1 

sup E[\ y f n \ 2 ] + e(£ i*ri a 5] + -c E(«r) 2 + -7 E(^ ,n ) 2 ) ^ ^, g , Ln , un . 

' i=o p0 i=o p0 i=o 

We reconsider Ito formula for |?/J' n | 2 , the take sup.,- before expectation. Using Burkholder- 

Davis-Gundy inequality for martingale part Yli=o ^'"^'"^^i+i' w ^h similar techniques, 
we get 

n-1 

E[snp \yf n \ ] < C e , ffl i« )£ /n + C„E[^2 \y P i ,n \ 2 S\ < C e , 9 ,L«,c/™ + C^Tsup E[\y?' n \ 2 ). 



j 



8=0 



It follows the desired results. □ 



4 Reflected Algorithms and their convergence 
4.1 Reflected Schemes 

This type of numerical schemes is based on reflecting the solution y n between two barriers by 
a n and k n directly. In such a way the discrete solution y n really stays between two barriers L n 
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and U n . After discretization of time interval, our discrete reflected BSDE with two barriers 
on small interval [tj, for < j < n — 1, is 

y- = y? + i + g(t jt %)6 + a] - k] - zjV~5e] +1 , (13) 

with terminal condition y™ = £ n , and constraint and discrete integral conditions hold: 

a] > 0, k?>0, cq-k? = 0, (14) 
L] < y?<U?, ( y J-L])a] = (yJ-U?)kJ = 0. 

Note that, all terms in (fTBl are ^"-measurable except y™ +1 and 

The key point of our numerical schemes is how to solve (y™, zj, a™, k™) from (fTBl using the 
(/™ +1 -measurable random variable y™ +1 obtained in the preceding step. First zj is obtained 
by 

Then ([THD with (HID becomes 



y? = £[yj +1 |^ n ] + g(t v yl z])8 + a] - k%, a] > 0, fc? > 0, (15) 
V] < y] <U?, ( y J-L])a] = (yJ-U?)kJ = 0. 

Set 6(2/) := 2/ - ^ y, 4)5. In view of (9(y) - O(y'), y - y'> > (1 - fy) |y - yf > 0, for 5 
small enough, we get that in such case 0(y) is strictly increasing in y. So 



y > L?^e(y)>e(i-), 

Then implicit reflected scheme gives the results with E[y^ +1 \QJ] = i(y™ +1 | e n =1 + y" + i| e r>=-i) 
as follows 

y] = e-\E[ y - +1 \g-] + a]-k]), 

k« = (E[y] +1 \G^+g{t h Ulz])5-U?) + , 

on the set {L™ < C/j 1 }, then we know that {yj — L™ = 0} and {y™ — C/j 1 = 0} are disjoint. 
So with (y? - L])a] = [y] - Uf)k] = 0, we have a] ■ k] = 0. On the set {L™ = Uf}, we 
get a? = (/?) + and k] = (/?)- by deflation, where I] := E\y* +X \G*\ + g(t j ,L],z^)5 - L]. 
So automatically a™ ■ k™ = 0. 

Our explicit reflected scheme is introduced by replacing y™ in (7 by i?[y™ +1 |{?™] in (fTBT) . So 
we get the following equation, 

y] = Eiy^gjj+git^Eiy^ig^z^S + a^-V;, a™ > 0, ^ > 0, (16) 
< V?<U?, (yJ-L])a^ = (yJ-U^ = 0. 
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Then with E[fJ +1 \Gf 



we get the solution 



V" 



T" 



(17) 



4.2 Convergence of Reflected Implicit Schemes 

Now we study the convergence of Reflected Schemes. For implicit reflected scheme, we 
denote 



[t/S] 



Y t —y[t/s\i z * — z \t/s\i a — / j a i ) A t 

and for explicit reflected scheme 



i=0 



[t/a] 



[t/a] 

i=0 



[t/a] 



n — n 



i — V[t/S\ ' Z t — z [t/<5] > 



i=0 



i=0 



First we prove an estimation result for (?y n , z n , a n , k n ). 

Lemma 4.1 For 5 such that 6(1 + 2/j, + 2// 2 ) < 1, there exists a constant c depending only 
on /i and T such that 



n-l 



F[sup|^|+^|^| 2 5 + 

J 3=0 



n-l 
j=0 



n-l 
3=0 



n JJn . 



Proof. First we consider the estimation of a™ and kf. In view of U\ < Y n n < Uf, we 

l l j j j 

have that 



a] < (E[L] +1 \g?) + g(tj, L], z^)6 - Lfj = 6(1, + gfa L], z]))' , 
kj < {E[Uf +1 \gj]+g(t v U;\z-)5-U-) + = 5{u J + g(t J ,U^z^y 

We consider following discrete BSDEs with = y™ = £ n , 



v? = y^ +1 + [g(t jl y^,^) + (l J )- + g(t v L],^)-}5-z^V~6e] +1 , 
Vl = y^ +1 + ig(t J ,y?^)-(u J )+-g(t J ,U^^)+]5-^V5e] +1 . 

Thanks to discrete comparison theorem in [20], we have y™ < y™ < yj, so 

£[sup \y™\ 2 } < max{£[sup \y"\ 2 }, E[sup \y"\ 2 }} < c C^ t9iL n jUn . 



(19) 



13 



The last inequality follows from estimations of discrete solution of classical BSDE (y™) 2 and 
(?7™) 2 , which is obtained by Ito formulae and the discrete Gronwall inequality in Lemma [XT! 
For Zj, we use 'discrete Ito formula' (cf. [20]) again for (y™) 2 , and get 

n— 1 n— 1 n— 1 n—1 

|2 



*=j *=J i=j i=j 

n— 1 n—1 -. n— 1 

< ^tin 2 + E °> °)i 2 6 + ^ + 2 ^ + E i^ n i 2 + 2 E i^ n i 2 ^ 

1 n— 1 n— 1 

+a£[sup((L?)+) 2 + sup((^)-) 2 ] + -E[(J2 S") 2 + (E k W> 



using (y™ — L")a™ = and (y™ — Uf)k™ = 0. And from (JIB"]) , we have 

n—1 n— 1 

E(E«") 2 < 4^E[^) 2 + ^(^,0,0) 2 + /i|^| 2 + ^|^| 2 ], 

n—1 n— 1 

i?(E fc i) 2 ^ 4^E[M 2 + ^0,0) 2 + ^|^| 2 + /i|^| 2 ]. 

J=i jr'=i 

Set a = 32/i, it follows 

n— 1 1 n—1 n—1 

%if + I E l^l 2 5 ] ^ ^[l^l 2 + C 1 + g^) E I^C**, 0. 0)| 2 5} + 5(1 + 2^ + 2/i 2 ) £ M 

i=j ^ i=j t=j 

n-1 



+32^£[sup((L?) + ) 2 + sup((^)-) 2 ] + ^E-£[(l 

3 

+ ^E[| L "| 2 +I C/ : 
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H 1 
j? I 



With ffT9l . we obtain Yli=o \ z i\ ^ — c ^^ n ,g,L n ,u n - Then applying these estimations to ([TBI , 
we obtain desired results. □ 

With arguments similar to those precede Proposition 13.11 the laws of the solutions 
(Y,Z,A,K) and (Y n , Z n , A n , K n ) or (Y n , Z n ,T , K n ) to reflected BSDE depend only on 
(P B ,r- 1 (P B ),<7, ■*I 1 (Pb),% 1 (Pb)) and (P B n, r^P^),,?, ^^P^)) where 

/ _1 (Pb) (resp. /"HPb")) is the law of /(B) (resp. /(B n )) for f = T, V 1} * 2 . So if we 
concern the convergence in law, we can consider these equations on any probability space. 

From Donsker's theorem and Skorokhod representation theorem, there exists a probabil- 
ity space satisfying sup 0<t<T \B™ — B t \ — >• 0, as n — > oo, in L 2 (JF T ), since Sk is in L 2+(S . And 
it is sufficient for us to prove convergence results in this probability space. Our convergence 
result for the implicit reflected scheme is as follows: 
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Theorem 4.1 Under Assumption \2.3\ and suppose moreover that g satisfies Lipschitz con- 
dition (QP, we have when n — > +00, 

E[sup\Y t n -Y t \ 2 } +E [ \Z? - Z t \ 2 dt ^ 0, (20) 
t Jo 

and A n t - K? -> A t -K t in L 2 (j^) ; forO<t<T. 
Proof. The proof is done in three steps. 

In the first step, we consider the difference between discrete solutions of reflect implicit 
scheme and of penalization implicit scheme introduce in section 4.1 and section 3.1, respec- 
tively. More precisely, we will prove that for each p, 

n—l 1 

£[sup \y] - yff] + SE^lz? - zf n \ 2 < cC^ g ^ Un —. (21) 
i j=o VP 

Here c only depends on /i and T. From ([9]) and ( |T3|) . applying 'discrete Ito formula' (cf. 

TO to (y?-yf n ) 2 , we get 

n-1 

p,n 1 2 



E\y?-y?' n \ 2 + 5Ej2\z? 

i=j 

n—l 

i=j 

n—l n—l 

+2eJ2m - yf n m - - 2eJ2m - yrm - Kn\ 



1=] i=j 

From (HU), we have 



(y? - vHW - O = ivi - l?K - {yT - l?x - (% n - W + (»T - ^K' n 

< (yr-L7r<-((yf' n -L-)-f, 

Similarly we have (yf — yf n ){k™ — kf ,n ) > —{Vi' n — U^)k^. By ( ITS]) and the Lipschitz property 
of g, it follows 

r n—l 




E\yJ-yf n \ 2 +°-Ej2\z?-zr\ 2 

i=j 



n—l n—l 

< (2a* + 2» 2 )8Ej2u - yn 2 } + 2Ej2M' n - m~< + (yr - wrm 

i=j i=j 

1 

n—l / n—l \ 2 / n—l 



< {2^ + 2^ 2 )5EY,M ~yn 2 ] + 2 ^£((2/T-^n 2 ^^((^+^,L«,^))- 



i=j \ i=j / \ i=j 
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1 1 
n—1 \ 2 / n— 1 \ 2 

i=j J \ i=j 

1 1 
n—1 rs / -i n—1 \ 2 / n—1 \ 2 



-^2 



/ n-l \ 2 / n-1 

Then by estimation results in Lemma 13.21 Lemma 14.11 and discrete Gronwall inequality in 
Lemma [3. 1[ we get 

n—1 1 

sup E \y] - yf f + 6E £ K - afT < c . 

i 8=0 

Apply B-D-G inequality, we obtain ( l2~li) . 

In the second step, we want to prove ( 120|) . We have 

£[sup \Y t n - Y t \ 2 ] +E[f \Z? - Z t \ 2 dt] 
t Jo 

< 3£[sup|lf -Y t \ 2 + [ \Z P -Z t \ 2 dt] + 3£[sup |>7 - FH 2 + / \Z? - Z?' n \ 2 dt] 

t Jo t Jo 

+3E[sup\Y t p -Y t p ' n \ 2 + [ \Z p t - Zf' n \ 2 dt] 
t Jo 

< 3Cp^ + cC^, g . L ^ Un p- L 2 +3E[snp\Y p -Y p ' n \ 2 + / \Z P - Z p ' n \ 2 dt], 

t Jo 

in view of ( 12T1) and Theorem 12.11 For fixed p > 0, by convergence results of numerical 
algorithms for BSDE, (Theorem 12 in |5] and Theorem 2 in [25]), we know that the last two 
terms converge to 0, as 5 — > 0. And when 5 is small dominated by £ n , 

g, L and U. This implies that we can choose suitable 5 such that the right hand side is as 
small as we want, so (12"0|) follows. 

In the last step, we consider the convergence of (A n , K n ). Recall that for < t < T, 

AZ-K" = Y Q n -Y t n - f g(s,Y s n ,Z:)ds+ f Z?dB?, 

Jo Jo 

[ g(s,Y p > n ,Z p > n )ds+ [ Z p s ' n dB n s . 
Jo Jo 



y[P> n _ j£Pi n — Y p,n — Y p ' n 



By (l2Til and Lipschitz condition of g, we get 



E[\(A? - K?) - (Ar - Knf] < c Cm 



1 



n JJn 
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Since 

E[\(A» - K?) - (A t - K t )\ 2 ] < 3E[\(A? - K?) - (A? n - Kf' n )\ 2 ] + 3E[\(A p t - K?) - (A t - K t 

+3E[\(A p t - Kf) - (A p t ' n - Kn\ 2 ] 

< c(c^ tU n + c^u)-^ + mui - Kf) - (Ar - Km 2 ] 

with similar techniques, we obtain — K™) — (A t — K t )\ 2 ] — > 0. Here the fact that 

(At' n — Kf ,n ) converges to (Af — Kf) for fixed p follows from the convergence results of 

(yr,zn to (y?,z?). □ 

4.3 Convergence of Reflected Explicit Scheme 

Then we study the convergence of explicit reflected scheme. Before going further, we need 
an estimation result for (y n ,z n ,a n ,k ). 

Lemma 4.2 For 5 such that <J(| + 2yU + 4/i 2 ) < 1, there exists a constant c depending only 
on fi and T , such that 



n-l 



E[sup\yf] + E[J2\z]\S + 



n-l 



3=0 



-rn 

3 



n—1 
3=0 



< C C^ri g Ln JJn 



Proof. We recall the explicit reflected scheme, which is: 



y 



L] < !,]< I-;. L?)a? = = 0. 

Then we have 



(22) 



\y 



n|2 



irin 1 2 I -n 1 2 



\y 



j+ii 



\Z.; 



6 + 2y] +l ■ g(tj, E[y] +1 \g-],z])5 + 2yJ-c%- 2y? • % (23) 



"1 ^«\|2r2 i-rfliVl 



+\g(t J ,E[y] +1 \g-],z])\ 2 s 



JJ + L\-'] J' ~J 

+2g(t j ,E[y] +1 \g^z])zyV5e] +1 - 2(a» -T%)%VSe 



3 ■ 
% 



(^) 2 - 2yp"V5e] +1 



->3 3 



_n 

■-j+v 



In view of (y™ — L™ )a™ = (y™ — U™)7c* = 0, a™ and fc™ > 0, and taking expectation, we have 



J > 3 



E\f;\ 2 + E\zJ\ 2 5 < E\y r ; +l \ 2 + 2E[y] +l ■ g(tj, E\fi +1 \g?],z$)]6 + 2E[(L n J ) + ■ a]] + E[{U]T ■ k 

+E[\g(t j ,E[y* +1 \g?],z»)\ 2 6 2 } 

A 



< m + i\ 2 + ( 5 + ^ 2 ) E l\9(t3, o, o)| 2 ] + qs + ^H 2 )E[{z^r 

+6(1 + 2/i + V + 3pL 2 S)E\y] +1 \ 2 + 2E[(L]) + ■ a 1 *] + E[(Uy)~ ■ k 
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Taking the sum for j = i, ■ ■ ■ , n — 1 yields 



1 n— 1 

e\v?\ 2 + ^e\z^s 



(24) 



n-l 



n-1 



< 



£|ff + {5 + 35 2 )£ £[|<?(t, ; 0, 0)| 2 ] + 5(1 + 2^ + V + 3fi 2 6)E £ |j£ 



j=i 



3=i 



1 n— 1 n— 1 

+a£[sup((L?)+) 2 + sup((f/;) + ) 2 ] + -E[(£^T + (£^) 2 ], 

i i ol . , 

3=l 3=l 

where a is a constant to be decided later. From (fTTj) . we have 

a] < {E[L] +1 \gj} + g( tj , Eifi +1 \g?],z?)8 - U])' = (/, + g( tj , E[yJ +1 \g-},zJ))S, 
V; < {E[U? +1 \g}] + g(t v E\fi +1 \g?],%)6 - U«) + = ( Uj + g(t 3 , E[y? +1 \g?\, z]))~5. 
Then we get 



n-l 



n-l 



j=i 

n-l 

E(£k 



a?) 2 < 



45E J2m 2 + git,, 0, 0) 2 + /* a (^ +1 |^]) a + A*] 



(25) 



n-l 



< 



4^ £[(^.) 2 + g(t 3 , 0, 0) 2 + /A^il^]) 2 + f(z? 



Set a = 32/i 2 in it follows 



n-l 



< 



1 " x 

j=i 

E\C\ 2 + ( 5 +^2 S + ^ 2 )E 'Y^l\g{t v 0, 0)| 2 ] + 32 yU 2 £[sup((L?) + ) 2 + sup((£/ : 



n-l 



7) + ) 2 ] 



n-l 



n-l 



+<K^ + 2/x + 4^ 2 + 3f6)E \V] + i\ 2 + ^E^U) 2 + «) 2 ]. 

Notice that 3/i 2 <5 < 1, so 3fi 2 5 2 < 5. Then by applying the discrete Gronwall inequality in 
Lemma 13. 1[ and the estimation of and kj follows from (125]) , we get 

n-l 



n-l 



snpE[\y^\ 2 ] + E[Y^\z r ;\ 2 5 + 

j 3=0 



3=0 



—n 

3 



+ 



n-l 

£ 



< C C^n g J^n^ 



We reconsider (123]) . as before take sum and sup^-, then take expectation, using Burkholder- 
Davis-Gundy inequality for martingale part, with similar techniques, we get 

n-l 

£[sup \y]\ 2 ] < Qn g , L «^ + C^ElY^ \y]\ 2 S< £[sup \y]\ 2 } < C^^u* + C^T sup E[\y]\ 2 ], 



3=0 
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which implies final result. □ 

Then our convergence result for the explicit reflected scheme is 



Theorem 4.2 Under the same assumptions as in Theorem 4-1. when n — > +00, 

£[sup|T7 -Y t \ 2 } + E [ \T t - Z t \ 2 dt -> 0. (26) 
t Jo 

And T t -K n t -> A - #t in L 2 (^i), forO<t<T. 

Proof. Thanks to Theorem 14.11 it is sufficient to prove that asn-> +00, 

n-l 

E[snp fiff - y'f] + E J2W ~ z?\ 2 8 ^ 0. (27) 
J i=o 

Since 

iff = y? +1 +g(t j: y^z?)8 + a]-k?-z?V8e] +1 , (28) 

we get 

= E \vj + i - ^ + i| 2 - §E Vl - ^l 2 + ^I(y" - vlMtj> yJ> Z D - g(ts, m^ll^M 

-E[5{g{t h ylz])-g{t h E[r j+l \Q]]^)) + (a] -a]) - (k] -^)] 2 
+2E[(y] - y])(a] - a])] - 2E[{$ - yj)(k? ~ 

< E \v? + i - y] + i\ 2 - SE \ z l -*]\ 2 + 2SE l(y] - PjMtj, y], 4) - gfa E \y n 0+l \Q-}X 3 ))} 

in view of 

(y?-y])(a]-a%) = [y] - L])a] + (y] - L])(<%) 

-(jff - L])a] - (jff - 
< 

[y] -T 3 ){k] -V]) = [y] - U»)q + (y] - UJ)V; 

> o. 



We take sum over j from i to n — 1, with £™ — £™ = 0, then we get 

n— 1 n— 1 



n—1 r- n— 1 





n _ -n|2 

z i z j I 



n-l 
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Since y] - E[yf +1 \Gf] = gfo, E[y? +1 \g?],z?)5 + a]- k], we have 



2ii8E[\y?-y]\-\y?-E[y] +1 \g 

2it6E[\y? -y]\- y] -y] + g(tj, Effi +1 \g*}), zj)5 + a^-k 



—n 

3 



< (2/x + 1)8E[\$ -y]\ ] + itSEpPQgitj, 0, 0)\ z + f \y] +1 \ + H> |^| + W + (* 
Then by Lemma [4.21 we obtain 



E \yJ-y!\ 2 + oY, E 



n-l 



z n - z n \ 

Z 3 Z 3 I 



< (2/x 2 + 2(1 + 1)6^2 E\y 3 l -r 3 \ +SC^ Ln>un . (29) 



j=% j=i 
By the discrete Gronwall inequality in Lemma 13.11 we get 

£ '/', - f; |J 

■vn-1 



supply? -y?| 2 < C5V 2 ^ 2+1 ) T . 



With (|2S|), it follows E[SJ2]Zo E \ z j ~ Zj\ ] < C^ 2 - Then we reconsider (J2SD, this time we 
take expectation after taking square, sum and sup over j. Using Burkholder-Davis-Gundy 
inequality for martingale parts and similar tachniques, it follows that 



n-l 



£sup \y? - y]\ 2 < CE^ E\yJ 5} < CT sup E \y] - T 3 



3< n 



j=0 



which implies (J2j 

For the convergence of (A , K ) , we consider 



K - K 



1 1 t 



g(s,Y n s X)ds+ / Zjdff 



AZ-K? = Y n -Y t n - f g(s,Y s n ,Z:)ds+ T ZyBl\ 

Jo Jo 

then the convergence results follow easily from the convergence of A n , ( |26l) and the Lipschitz 
condition of g. □ 



5 Simulations of Reflected BSDE with two barriers 

For computational convenience, we consider the case when T — 1. The calculation begins 
from y™ = £ n and proceeds backward to solve (y™, Zj, a™, fc"), for j = n—1, •••1,0. Due to the 
amount of computation, we consider a very simple case: £ = $(£>i), L t = ipi(t, B(t)), U t = 
?p2(t, B(t)), where $, ipi and ip2 are real analytic functions defined on R and [0, 1] x R respec- 
tively. As mentioned in the introduction, we have developed a Matlab toolbox for calculating 
and simulating solutions of reflected BSDEs with two barriers which has a well-designed inter- 



face. This toolbox can be downloaded from http:/ /159.226.47.50:8080/iam/xumingyu/English/index.jsp 
with clicking 'Preprint' on the left side. 
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We take the following example: g(y, z) = — 5 \y + z\ — 1, $(x) = \x\, ^i(t,x) = — 3(x — 
2) 2 + 3, * 2 (t, x) = (x + l) 2 + 3(t - 1), and n = 400. 

In Figure 1, we can see both the global situation of the solution surface of y n and its 
partial situation i.e. trajectory. In the upper portion of Figure 1, it is in 3-dimensional. The 
lower surface shows the barrier L, as well the upper one is for the barrier U. The solution 
y n is in the middle of them. Then we generate one trajectory of the discrete Brownian 
motion (-B™)o<j<n ; which is drawn on the horizontal plane. The value of y™ with respect to 
this Brownian sample is showed on the solution surface. The remainder of the figure shows 
respectively the trajectory of the force A 7 - = J2i=o ^ an d Kj = Yli=o k?- 




t t 
Figure 1: A solution surface of reflected BSDE with two barriers 

The lower graphs shows clearly that A n (respective K n ) acts only if y n touches the lower 
barrier L n , i.e. on the set {y n = L n } (respective the upper barrier U, i.e. on the set 
{y n = U n } ), and they never act at the same time. 

In the upper portion we can see that there is an area, named Area I, (resp. Area II) 
where the solution surface and the lower barrier surface (resp. the solution surface and the 
upper barrier surface) stick together. When the trajectory of solution y™ goes into Area I 
(resp. Area II), the force A™ (resp. K^) will push y™ upward (resp. downward). Indeed, if 
we don't have the barriers here, yj intends going up or down to cross the reflecting barrier 
L™ and [7™, so to keep staying between L™ and [7™, the action of forces A™ and are 
necessary. In Figure 1, the increasing process A™ keeps zero, while increases from the 
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beginning. Correspondingly in the beginning y" goes into Area II, but always stay out of 
Area I. Since Area I and Area II are totally disjoint, so A™ and K™ never increase at same 
time. 

About this point, let us have a look at Figure 2. This figure shows a group of 3- 
dimensional dynamic trajectories (t,-, B™, YJ 1 ) and (tj, BJ, ZJ), simultaneously, of 2-dimensional 
trajectories of (tj, YJ 1 ) and (tj, ZJ). For the other sub-figures, the upper- right one is for the 
trajectories AJ, and while the lower-left one is for KJ", then comparing these two sub-figures, 
as in Figure 1, {aj ^ 0} and {kj ^ 0} are disjoint, but the converse is not true. 

trajectory of y(t) in 3-D trajectory of z(t) in 3-D 




0.5 1 0.5 1 0.5 

t t t 

Figure 2: The trajectories of solutions of 



Now we present some numerical results using the explicit reflected scheme and the 
implicit-explicit penalization scheme, respectively, with different discretization. Consider 
the parameters: g(y, z) = — 5 \y + z\ — 1, = ^i(t, x) = — 3(x — 2) 2 + 3, x ) — 
(x + l) 2 + 3t-2.5: 

n = 400, reflected explicit scheme: = —1.7312 



penalization scheme: 
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n = 1000, reflected explicit scheme: y$ = —1.7142 



penalization scheme: 
n = 2000, reflected explicit scheme: 

penalization scheme: 
n = 4000, reflected explicit scheme: 

penalization scheme: 
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From this form, first we can see that as the penalization parameter p increases, the 
penalization solution y^' 71 tends increasingly to the reflected solution y^. Second, as the 
discretaization parameter n increases, the differences of y$ with different n become smaller 
as well as that of y p,n . 



6 Appendix: The proof of Theorem 12.1 



To complete the paper, here we give the proof of Theorem 12 . 11 

Proof of Theorem 12.11 (a) is the main result in [16] . So we omit its proof. 

Now we consider (b). The convergence of (Y p , Z p ) is a direct consequence of [23]. For 
the convergence speed, the proof is a combination of results in [TH] and [23] . From |16j . we 
know that for ©, as m -> oo, Y t m ' p / Y P in S 2 (0,T), Z™' p -> Z» in l£(0,T), A™' p -> A P 
in S 2 (0,T), where (Y P , Zf, A P ) is a solution of the following reflected BSDE with one lower 
barrier L 

-dY P = 9 (t,n,Z P )dt + dA P -p(Y p -U t ) + dt-Z P dB t ,Y P n = ^ (30) 

Yf > L u [\y p - L t )dM = 0. 
J o 

Set K P = JqP(Y p - U s ) + ds. Then letting p -> oo, it follows that Y_ P \ Y t in S 2 (0,T), 
Z% — > Z t in L|r(0, T). By comparison theorem, gL4£ is increasing in p. So /* At, and 
< sup 4 [A P+1 - £ t ] < £ T +l - A T . It follows that A P -> A t in S 2 (0, T). Then with Lipschitz 
condition of g and convergence results, we get Iff — > in S 2 (0, T). Moreover from Lemma 
4 in [16J, we know that there exists a constant C depending on £, 0, 0), fi, L and £7, such 
that 



£[ sup in - yi 2 + I m- z t \ 2 )dt < 

0<t<T Jo y/P 

Similarly for ©, first letting p -> oo, we get Y t m ' p \ in S 2 (0, T), -> in L^(0, T), 
K rn, P ^ -^m . n S 2 (0)T); where (Y^^^K t ) is a solution of the following reflected BSDE 
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with one upper barrier U 

-dY? = g(t,Y^X)dt + m(L t -YT) + dt-dX l -ZTdB t ,Y^ = ^ (31) 

FT < U t , I (FJ - U t )dK? = 0. 
./o 



In the same way, as m — > oo, Y t /* Y t in S 2 (0, T), Z t — > Z t in L^-(0, T), and (v4 t , if t ) — > 
(A tl K t ) in (S 2 (0,T)) 2 , where = j^m(L s — Y^) + ds. Also there exists a constant C 
depending on £, 0, 0), /i, L and {7, such that 



E[sup \Y™-Y t \ 2 + / (ir" - Z t \ 2 )dt < 

0<t<T Jo 



c 



m 

T7P 



Applying comparison theorem to (jSJ) and fl30|) . (jUJ) and fl3Tj) (let m = p), we have F t < 
Yf < Yf. Then we have 

£[sup \Yf-Y t \ 2 ]<^, 

0<t<T y/P 

for some constant C. To get the estimate results for Z p , we apply Ito formula to \Yf — Y t \ 2 , 
and get 

E\Y Q P -Y Q \ 2 + \e f \Z P S -Z s \ 2 ds 
* Jo 

= fa + 2fi 2 )E [ \Y S P - Y s \ 2 ds + 2E [ (Yf - Y s )dA p - 2E [ (Y p - Y s )dA s 
Jo Jo Jo 

-2E [ (Y*-Y 8 )dK* + 2E [ {Yf — Y s )dK s . 
Jo Jo 

Since 



2E [ (Y s p - Y s )dA p = 2E [ (Y p - L s )dA p -2E f (Y s - L s )dA : 
Jo Jo Jo 

< 2pE ! (Y p -L s )(Y p -L s yds<0 
Jo 

and 2E J^(Y P - Y s )dK p > 2pE J*(Y P - U S )(Y P - U s ) + ds > 0, we have 

E / \Z p - Z s \ 2 ds < —, 
Jo Vp 

in view of the estimation of A and K and the convergence of Y p . 
Now we consider the convergence of A p and K p . Since 



A t -K t = Y — Y t — / g(s,Y s ,Z s )ds+ / Z s dB s , 

Jo Jo 

Al-Kf = Y p -Yf- f g(s,Y p ,Z p )ds+ f ' Z p dB s 

Jo Jo 
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from the Lipschitz condition of g and the convergence results of Y p and Z p , we have imme- 
diately 



E[ sup [(A t - K t ) - (A p - K p )] 2 ] 

0<t<T 

< BE[ sup \Y t p - Y t \ 2 + Afi [ \Yf -Y s \ 2 ds + C [ \Z P S - Z s \ 2 ds] < 

0<t<T Jo Jo \/V 

Meanwhile we know E[(A P T ) 2 + (K^) 2 } < oo, so A p and K p admits weak limit A and K in 
S 2 (0, T) respectively. By the comparison results of Y%, Y t p and Y_ P , we get 

dA p = p(Y p -L t )-dt<p(Y P -L t )-dt = d~A P , 
dK p = p(Y p -U t ) + dt>p(Y p -U t ) + dt = dK p . 

So dA t < dA t and dK^> dKj, it follows that dA t — dK t < dA t — dK t ^ On the other hand, the 
limit of Y p is Y, so dA t — dK t = dA t — dK t . Then there must be dA t = dA t and dK t = dK t , 
which implies A t = A t and K t — K t . □ 
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